Towards a bulk theory of fiexoelectricity 
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Fiexoelectricity is the linear response of polarization to a strain gradient. Here we address the 
simplest class of dielectrics, namely elemental cubic crystals, and we prove that therein there is no 
extrinsic (i.e. surface) contribution to fiexoelectricity in the thermodynamic limit. The flexoelectric 
tensor is expressed as a bulk response of the solid, manifestly independent of surface configurations. 
Furthermore, we prove that the flexoelectric responses induced by a long-wavelength phonon and 
by a uniform strain gradient are identical. 
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Fiexoelectricity is by definition the linear response of 
polarization P to strain gradient; it is measured by the 
fourth-rank Cartesian tensor n a pjs, i-e. 

Pa = Val3~<S , (1) 

where summation over repeated indices is understood; 
such tensor property is symmetry allowed in any 
dielectric. In recent years, there has been much interest 
in creating piezoelectric composites from materials 
which are not themselves piezoelectric, by exploiting 
fiexoelectricity [1]. There have been also some 
measurements [2], atomistic calculations [3], and even 
(very recently) first-principle calculations for BaTiC>3 
and SrTi0 3 [4]. However, the basic issue whether 
fiexoelectricity is a pure bulk effect — or instead it has 
a contribution which depends on the surface conditions 
of the sample — is unsettled to date, even for the simplest 
of crystalline dielectrics. It is worth mentioning that the 
analogous issue about the simpler case of piezoelectricity 
is far for being trivial either [5], and spurred much 
discussion among theorists until two decades ago [6]. 
The issue whether a material property is a bulk effect 
is of the utmost importance for theorists, since in 
the affirmative it can be in principle computed within 
periodic boundary conditions: for piezoelectricity, this is 
happening routinely since 1989 [7]. The most quoted 
theoretical paper on fiexoelectricity, by Tagantsev [8], 
was published in 1986. Tagantsev's message is opposite 
to the present one, even in the abstract; a brief critique 
of Tagantsev's work is provided at the end of this work. 

The present work is a first step towards a bulk theory 
of fiexoelectricity. Here we limit ourselves to the simplest 
possible case: elemental cubic crystals, having a primitive 
lattice: we prove that fiexoelectricity is indeed a genuine 
bulk property therein: there is no surface contribution to 
the flexoelectric response. The class of crystals addressed 
here includes solid rare gases and little more; therefore 
the present result looks having only academic interest. 
However, we believe that the additional effects (due to 



sublattice displacements and ionic screening) which occur 
in nonprimitive lattices can be also tamed at the price of a 
much clumsier algebra. This is briefly discussed towards 
the end of this work. 

The definition of Eq. (1) is incomplete, since one has to 
specify the macroscopic field E. It is customary to define 
the material constants (such as Born effective charges, 
piezoelectric and pyroelectric constants) as derivatives of 
P with respect to the relevant variable, taken at zero 
field. If we make such a choice in Eq. (1), then the 
flexoelectric polarization in a nonzero field E is then 

de^s £oo — 1 „ ,r,s 

Orp 47T 

where we have exploited cubic symmetry and is 
the dielectric constant (sq — £00 i n primitive lattices). 
Whenever /j, is a bulk effect, then Eq. (2) does not 
depend on the shape of the sample; the only effect of 
the shape is represented by a depolarization field and 
included in E [9, 10]. For a given shape (free-standing 
in zero external field), the depolarization field dictates 
the surface charge at the boundary, via Gauss theorem. 
Conversely, when a genuine surface effect contributes 
to the apparent polarization, some extra charge must 
accumulate at the boundary. Here we choose the worse 
possible shape (see below) , where the depolarization field 
is maximum, and we show that even in this case there is 
no extra boundary charge. The absence of surface effects 
is confirmed by the case (discussed below) of a crystal 
with no boundary and a built-in long- wavelength phonon. 

In a crystal of cubic symmetry the tensor fiap^s 
depends on three independent elements: these can 
be evaluated in principle by considering the uniaxial 
response — where strain, strain gradient, and polarization 
are all aligned — along three different directions, e.g. 
(100), (110), and (111). For the sake of simplicity, in 
the following we label as "a;" any of these directions, and 
we address the uniaxial flexoelectric coefficients, i.e. 

p 9e xx — 1 . . 

*x l^xxxx ~ a 1 Z ^x- \3) 
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FIG. 1: Slabs with a built-in uniform strain gradient 
(uniaxial); the x-direction is indicated by the arrow in each 
panel. Top: "transverse" case. Strain, strain gradient, 
and polarization are parallel to the slab, while the field 
vanishes. Bottom: '"longitudinal" case. Strain, strain 
gradient, polarization, and depolarizing field are all normal 
to the slab. Similar figures can be drawn even when x is not 
along a cubic axis. 



We consider a sample in the form of a slab, where 
de xx /dx is either parallel or normal to the slab: this 
is schematically shown in Fig. 1. In the former case E 
vanishes and P is parallel to x ("transverse"); in the 
latter case both E and P are parallel to x, and are 
related by E = -4ttP ("longitudinal") [9, 10]. We 
stress that such properties are a consequence of cubic 
bulk symmetry, even if x is not a cubic direction (Fig. 1 
only shows cubic directions for the sake of simplicity). 
Eventually, we will address infinitely thick slabs: the two 
thermodynamic limits are different because of the long 
range of Coulomb interaction. 

We focus on the longitudinal case: replacement of 
E = -4ttP into Eq. (3) yields 



— — 



9e, 



dx 



(4) 



We are going to use this key relationship in reverse, i.e. 
we will prove that: in our longitudinal geometry the 
macroscopic field (i) does not depend on what happens 
at the slab surfaces; (ii) is constant in the bulk of the 
slab, and (iii) is a well defined linear response function 
of the material. This amounts to prove that fi xxxx is a 
bulk material property, and that no surface contribution 
enters Eq. (1) for the class of solids dealt with here. 

Our key ingredient is the microscopic electric field E(r) 
linearly induced by a unit displacement of the £-th 
nucleus in the otherwise unperturbed lattice. Because of 
translational invariance, 



dE p (r) 



£ a p{r- R e ) 



(5) 



it goes to zero as an inverse power of |r| and admits a 
multipolar expansion. We remind that the x-direction is 
not necessarily a cubic axis. The induced charge is given 
by V • E = Airp; in a primitive lattice it is an odd function 
of r, and its dipole is zero because of the acoustic sum 
rule [11]. Therefore the induced perturbation is octupolar 
to leading order. 

We focus on the atomic plane whose equilibrium x- 
coordinate is X m = md, where d is the interplanar 
distance and m is an integer, and we reduce our 
problem to an effective one-dimensional one by averaging 
everything in the yz plane. The averaged microscopic E 
field in the x-direction induced by a rigid displacement 
u m of the m-th atomic plane in the x direction is 



E(x) — u m £ xx {x X m ), 



(6) 



where £ xx {x) is the average of £ xx (r) in the planes normal 
to x. The yz- averaged induced charge is then, by Poisson 
equation 

The conditions of zero induced monopole, dipole, and 
quadrupole lead, after integration by parts, to 

/GO 
dx £ xx (x) = 0, 
-OO 

dx x£ xx (x) = 0, (8) 
while the second moment of £ is essentially the octupole: 

(9) 



J — ( 



/ 

J — c 



dx x £ xx (x) Q XX xx. 



where Qxxxx is the third moment of the induced charge. 
Actually, al I moments of £ xx (x) are convergent integrals, 
which proves that one-dimensional electrostatics is short 
range: £ xx (x) — > faster than any polynomial for x — > 

±00. 

If our slab is subject to a constant strain gradient along 
x, the displacement of the m-th plane can be written as 
u m = ridm 2 /2. In fact the local strain at plane m, defined 

as 



^xx (X m ) — 



Um+1 — Um-1 
2d 



(10) 



is thus e xx (X m ) = rjm (77 dimensionless constant). The 
induced microscopic field, averaged over yz, is after 
Eq. (6): 



r)d 
Y 



S ( x )='-W m 2 £ xx (x-X m ). (11) 



and £ a/ 3 (r) is a well defined linear response of the system; 



mGslab 

Here we have assumed that the terms with X m close 
to the slab edge are identical to those in the bulk; 
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actually, they are different, but the surface effects (due 
to a nonextensive set of m's) cannot propagate deep in 
the bulk. Owing to the short range of the perturbation 
within our chosen geometry, extrinsic (i.e. surface) effects 
are ruled out. Actually we can replace the finite sum with 
the infinite one, over all m from — oo to oo.: 



E(x) 



rjd 
~2 



OO 

E • 

m— — oo 



&xx {% X m ) . 



(12) 



We stress that the logics leading to our formalism for 
the infinite solid — via slabs of growing thicknesses — 
amounts to performing the thermodynamic limit in two 
steps: first along yz and then along x. This guarantees 
that whenever the u m 's in Eq. (6) lead to a constant 
macroscopic field E, the corresponding macroscopic 
polarization obeys E = — 47rP. 

So far we have arrived at the microscopic field induced 
by a uniaxial flexoelectric distortion, Eq. (12), which is a 
wildly oscillating nonpcriodic function of x. At this point 
we need to address the macroscopic average of E(x), in 
order to prove that it is indeed a constant everywhere in 
the sample. According to textbooks [12] the macroscopic 
field is the convolution 



(E(x)) 



dx' w(x-x')E(x') = {w*E)(x), (13) 



where w(x) is a real function, nonzero in some 
neighborhood of x = 0, and normalized to unity; we are 
going to show that the function E(x) in Eq. (12) becomes 
a constant after convolution with a suitable w(x). We 
start with the standard window convolution function 



w 1 {x) = ±0(d/2-\x\ 



(14) 



where i? is the step function. Obviously w\* f extracts 
the average from any function f(x) periodic of period d. 
Next we define the family of convolution functions 

w 2 {x) = {w 1 * wi)(x), . . . w n (x) = (wi * w n _i)(x). (15) 

Suppose that f(x) is a polynomial of degree n times a 
periodic function (of period d) whose average is zero: it 
is easy to show that (w n+ i * f)(x) vanishes in this case. 
We write the microscopic field of Eq. (12) identically 



as 



E(x) 



' ^(a: — md) 2 £ xx {x — md) 

m 

( ^ (x — md) £ xx (x — md) 



2d 

■qx 



+ -tttI^ £xx(x-md). 



2d 



(16) 



If we now choose 1^3(2;), as defined above, to perform 
the macroscopic average, it is easy to see that both the 



second and third lines of Eq. (16) yield a vanishing result. 
The term in the first line is a periodic function, whose 
average is obviously constant 



Ey 



(E(x)) 



V 
2d 2 



f 

J — ( 



dxx 2 £ (x)-- 4n ^O^ 

(17) 

where Eq. (9) has been used. Since our dimensionless r) 
equals dde xx /dx, Eq. (4) yields 



£ O i3) 

6d : 



(18) 



which is the central result of this work, proving that the 
uniaxial flexoelectric coefficient can be expressed in terms 
of bulk linear response quantities. 

Using the same path as above, it is immediate to 
verify that constant strain — i.e. u m = e xx X m — induces a 
vanishing macroscopic field and a vanishing polarization, 
as it must be: a primitive lattice is nonpiezoelectric. 

In order to confirm our central result, we consider 
next a long-wavelength longitudinal phonon, generaliz- 
ing Martin's approach [5], originally devised for piezo- 
electricity, to the flexoelectric case. For a purely longitu- 
dinal phonon along x of amplitude u and wave vector k 
the displacements of the atomic planes are u m — u e lkXm , 
hence the induced microscopic electric field, after yz av- 
erage, is 



E(x) 



CO 

E 



ikX„ 



= uc ikx J2 e- ik ^- x ^£ xx (x-X m ).(19) 

m— — 00 

This is the product of a slowly varying envelope function 
times a periodic function. The macroscopic field and 
polarization at wavelength k are therefore 

00 

(E k ) = u(J2 e"^-*-) £ xx (x - X m )) 

Tfl — — OO 

(P k ) = = - 4^/°° dx e~ ikx £ xx (x), (20) 

Its lowest-order k expansion yields 

k 2 



(Pk) * 



Aird 2 



dx X &xx (*^) 



= _!^!n(3) 
6d ^ xxx 



(21) 



where Eq. (9) has been used; terms of order zero and one 
in k vanish owing to Eq. (8). At this point we remind 
that the macroscopic strain induced by a long- wavelength 
acoustic phonon is [5] e xx — iuk, hence its gradient is 
ike xx = —uk 2 . Eq. (21) becomes 



(Pk) 



0(3) 
^zz,xxxx 



de xx 



6d dx 



(22) 
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Comparison with our central result, Eq. (18), proves 
that the flexoelectric polarization induced by a long 
wavelength phonon and the one induced by a uniform 
strain gradient are the same. The trivial factor simply 
accounts for the fact that our primary \i definition is 
transverse and not longitudinal — see Eqs. (3) and (4). 

The present theory is strongly inspired by Martin's 
theory of piezoelectricity [5], although microscopic 
induced fields are addressed here instead of induced 
charges. In a piezoelectric crystal the leading multipolcs 
are dipoles and quadrupoles: according to Martin these 
uniquely determine the piezoelectric tensor. In the 
simple class of crystals considered here dipoles and 
quadrupoles are both zero: we have shown that octupolcs 
uniquely determine the flexoelectric tensor. In this sense, 
flexoelectricity can be regarded as the next higher order 
analogue of piezoelectricity. 

The above considerations also hint at how the present 
theory can be extended to real materials of interest, 
for instance to cubic perovskites [2, 4]. In such crystal 
structure our basic linear-response ingredients, Eq. (5), 
acquire a sublattice index, i.e. £ a p(r) — > £ SjCt p(r); the 
dipoles no longer vanish but their sum over s does, 
owing to the acoustic sum rule [11]; the quadrupoles 
vanish since each nuclear site is centrosymmetric. The 
above derivation can therefore be extended to this 
case, taking into account the fact that the dipoles are 
coupled to the internal strain and yield a nonvanishing 
contribution, similarly to what done by Martin in the 
case of piezoelectricity [5, 13]. In a nonprimitive lattice 
a related qualitative feature also occurs: Sq 7^ Sec. 
According to general macroscopic considerations [9, 10] 
the longitudinal and transverse polarizations differ by 
a factor Eq, not [14]. This implies that the 

perturbed nuclear coordinates in the longitudinal and 
transverse geometries are different, at variance with 
the elemental case sketched in Fig. 1; this difference 
amounts nonetheless to an internal strain, whose effect 
can be tamed (once more, in the same way as in 
piezoelectricity [13]). 

The present theory is also indebted to the theory of 
"absolute deformation potentials" [15], which addresses 
the electrostatic potential lineup between two differently 
strained regions of the same solid. Suppose we have a 
"homojunction" normal to x, i.e. a region on the left 
where the strain is zero and one on the right where the 
strain is e xx . The two regions are joined by a slab where 
the strain is graded; e xx coincides with the integrated 
strain gradient across the slab. In the transition slab 
the flexoelectric response is nonzero, and a depolarizing 
field is present (given the longitudinal geometry); its 
integrated value coincides with the potential lineup. Not 
surprisingly, the theory of Ref. [15] is based on induced 
octupoles as the present one is. 

Last but not least, we strongly disagree with the 
treatment of flexoelectricity provided by Tagantsev in 



Ref. [8], and subsequently adopted by other authors 
(e.g. Ref. [3]). This treatment is based on the rigid- 
ion model throughout; since a primitive lattice has 
zero ionic charges, the model is inadequate to address 
flexoelectricity in our case study. But there is more 
to say: electrons behave quite differently from classical 
point charges, and quantization phenomena dominate 
macroscopic polarization, sometimes in counterintuitive 
ways [16]. One of the purportedly non-bulk contributions 
to flexoelectricity is the charge second moment per unit 
volume of the unperturbed macroscopic sample (symbol 
"J" in Refs [3, 8]). For a centrosymmetric insulating 
material, whose surfaces are also insulating, this quantity 
vanishes owing to the theorem of quantization of the 
surface charge [17]; the theorem is generally violated 
by a system of classical point charges. Other flaws of 
Ref. [8] are the presence of conditionally convergent sums, 
without any prescription about their thermodynamic 
limit, and the lack of any relationship between shape and 
macroscopic field in the polarized sample. 

In conclusion, we have achieved a first basic step 
towards a bulk theory of flexoelectricity. We have 
shown that — in the simplest class of dielectrics at least — 
there is no surface contribution to flexoelectricity, which 
is a therefore a pure bulk effect. Furthermore the 
bulk flexoelectric responses for either long-wavelength 
phonons or uniform strain gradient are the same; they 
are expressed in terms of quantities which are manifestly 
surface-independent . 
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